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The penta-quark(5Q) baryon is studied in anisotropic quenched lattice QCD with renormalized 
anisotropy as/at=4 for a high-precision mass measurement. The standard Wilson action at /3 = 5.75 
and the 0{a) improved Wilson quark action with k=0. 1210(0.0010)0. 1240 are employed on a 12=^x96 
lattice. Contribution of excited states is suppressed by using a smeared source. We investigate both 
the positive- and negative-parity 5Q baryons with 7 = and spin J = 1/2 using a non-NK-type 
interpolating field. After chiral extrapolation, the lowest positive-parity state is found to have a 
mass, me = 2.25 GeV, which is much heavier than the experimentally observed 0^(1540). The 
lowest negative-parity 5Q appears at me = 1.75 GeV, which is near the s-wave NK threshold. To 
distinguish spatially-localized 5Q resonances from NK scattering states, we propose a new general 
method imposing a "Hybrid Boundary Condition (HBC)'\ where the NK threshold is artificially 
raised without affecting compact five-quark states. The study using the HBC method shows that 
the negative-parity state observed on the lattice is not a compact 5Q but an s-wave NK-scattering 
state. 

PACS numbers: 12.38.Gc, 12.39.Mk, 14.20.-c, 14.20.Jn 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) is the fundamental 
theory of hadron physics. It specifies elementary inter- 
actions at the level of quarks and gluons, which serve 
as building-blocks to construct complicated interactions 
among hadrons. At high energy, the QCD coupling con- 
stant diminishes due to asymptotic freedom, and pertur- 
bative QCD works well. At low energy, however, with 
the growth of the QCD couphng constant, perturbation 
theory breaks down, and the system is dominated by non- 
perturbative effects. 

One of the most important features of non- 
perturbative QCD is color confinement, i.e., all hadrons 
are formed as color singlet combinations of quarks and 
gluons. In the conventional picture, mesons consist 
of a quark and an antiquark, while baryons are made 
of three quarks. QCD, however, allows more general 
hadrons, often called exotic hadrons, which contain extra 
pairs of quark and antiquark, or gluons. Possible candi- 
dates include multiquark hadrons, glueballs, and hybrid 
hadrons. Experimental searches found no manifestly ex- 
otic hadrons before the year 2002. Therefore, the discov- 
ery of the manifestly exotic narrow resonance 8"'"(1540) 
by LEPS group at SPring-8 was an epoch making event 
in the hadron physics 1] . 

O"'' is confirmed to have baryon number B — I, charge 
Q = +1 and strangeness S = +1. It has to contain 
at least one anti-strange quark and its simplest config- 
uration is uudds. Therefore, O'^ is a manifestly exotic 
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penta-quark state. 

The experimental search fll] of 0+ was motivated by 
a theoretical prediction given in Ref. 2J. It is noted that 
0+ had been considered several times in the history 0, 
0, 0] . The discovered peak in the nif + invariant mass is 
centered at 1.54±0.01 GeV with a width smaller than 25 
MeV. A number of groups 0,0,0,01 confirmed the LEPS 
discovery, among which Ref. [3 claims that must be 
isoscalar, since no G+"'" is observed in the pK'^ invariant 
mass spectrum. 

Numerous theoretical studies of penta-quark baryons 
have appeare d since its discoverv l a [lo| . [ill ll^ 115 Il3 . 
Illfll WaS; 19, 20, 21, 22, 23, M,^5^M7'^' " 
l3ll l32l l33l l34. 35,. . One of the important issues is the 
spin and parity of O"*". Although J'^ = is suggested 
in the original work of Ref.|2| based on the chiral soliton 
model, there is no direct experimental information avail- 
able on spin and parity [33 • An educated guess for the 
spin is J = 1/2, since the color- magnetic interaction fa- 
vors smaller spin states in general. There is, however, no 
consensus on parity at all. Experimental determination 
of the parity of 0+ happens to be quite challenging |38| . 
while opinions are divided in theory );39, 4Qj. 

Positive parity is supported by various model calcu- 
lations, i.e., sohton models "5, 1 [H [H [UliJ , chiral 
bag model J^], the Jaffe-Wilczek's diquark model [l3 |. 
the Karliner-Lipkin's diquark-triquark model [T5l |. some 
quark model calculations E [Hill E 113 , and other 
model calculations ^21, ■ Negative parityis supported 
by some other quark model calculations 0, m llli lU; (2^ j 
and QCD sum rules |2||. 

There are several lattice QCD studies of penta-quarks 
available HI E H EE El 113, but they have not 
reached a consensus yet. Except for Ref. 01 j all these cal- 
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culations support that negative parity states are lighter 
than the positive parity ones and that positive-parity 
penta-quarks may appear at the mass above 2 GeV. 
Refs .143.1431 claimed existence of a bound negative-parity 
penta quark, whereas Ref.^^ has observed no evidence 
for narrow resonances. Naively, the negative result by 
Ref.'45'l may be understood as a natural consequence of 
their choice of the NK-type interpolating field, since it is 
expected to couple strongly to the NK continuum state. 
Therefore, under these circumstances, it is important to 
use various interpolating fields and accumulate more data 
to give a precise prediction. 

The aim of this paper is twofold, (i) We provide a 
high-precision data on the penta-quark baryon 0+ by 
using the anisotropic lattice QCD. (ii) We propose a new 
method, which enables us to distinguish compact penta- 
quarks from NK scattering states and apply it to the 
negative-parity 0"*". 

The paper is organized as follows. In Sect. ^ we dis- 
cuss our choice of the interpolating field for the penta- 
quark 0''", and introduce the parity projection. We 
adopt the standard Wilson gauge action at /3 — 5.75 
on the 12"^ x 96 lattice with the renormalized anisotropy 
^ = as /at = 4. The use of the anisotropic lattice is 
known to serve as a powerful tool for high-precision mea- 
surements of temporal correlators (48l l49l l5d . |5l| . For 
quark action, we adopt 0(a)-improved Wilson (clover) 
action with four values of the hopping parameter as 
K = 0.1210(0.0010)0.1240. We also employ a smeared 
operator to enhance the low-lying spectra. 

In Sect, mil we present our numerical results for both 
positive and negative parity channels. For both cases, we 
observe a rather stable plateau in the effective mass plot. 
For the positive parity case, we obtain, after the chiral 
extrapolation, me = 2.25 GeV, which is much heavier 
than the observed O+(1540). For the negative parity 
case, we obtain me = 1.75 GeV, which is rather close 
to the empirical value. Our data thus suggest that the 
negative-parity state is the ground state of the penta- 
quark spectrum. 

To clarify whether this negative parity state is really 
a localized resonance or not, we propose a new general 
method with a new ^^hybrid boundary condition (HBCf^ 
for the quark fields. Sect. lIVl is devoted to description 
of the HBC method. In the HBC method, the spatial 
boundary condition on the quark fields are artificially 
taken so that N and K must have non-vanishing spatial 
momenta, while it does not affect localized penta-quark 
states. As a consequence, the HBC raises the s-wave NK 
threshold by a few hundred MeV without changing the 
mass of a compact resonance state, 8+. 

By applying the HBC method, we investigate the 
negative-parity state further. Particularly, we exam- 
ine difference of the spectra between the two bound- 
ary conditions. We find that, in the effective mass 
plots, the plateau is raised by about the same amount 
as the shift of the NK-thresholds. We therefore con- 
clude that there is no localized resonance state below 
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the state observed under the standard boundary condi- 
tion turns out to be an NK scattering state. 

In Sect. a summary is given and possible implica- 
tions of our results are presented. 

II. INTERPOLATING FIELD AND PARITY 
PROJECTION 

We consider a non-NK type interpolating field for the 
penta-quark 0+ as 

Oa = ^abceadeebfg (u^Cj^de) (uJCdg) {Cs'^)a, (1) 

where a denotes the Dirac index, and a — g are color in- 
dices. C = 7472 denotes charge conjugation matrix. The 
quantum number of Eq. ^ is spin J = 1/2 and isospin 
/ = 0. It has an advantage that it cannot be decom- 
posed into a product of N and K in the non-relativistic 
limit. Hence, we may expect that its coupling to NK 
scattering states is rather weak. This is the crucial dif- 
ference between our calculation and that of Ref . (i^l , who 
adopted the NK type interpolating field, which is ex- 
pected to couple to NK scattering states strongly. The 
NK type interpolating field would not be an appropriate 
choice in studying the narrow penta-quark resonance in 
this sense. 

We present Fierz rearrangements of Eq. in Ap- 
pendix. ^ We can explicitly see that the NK compo- 
nent in the correlation function is suppressed by a factor 
1/16 compared to that from the NK-type operator. This 
non-NK type interpolating field Eq. was adopted in 
Refs. SIM Bl- 
under the spatial reflection, the quark fields transform 
as q{t,x) —> 74g(t, — x). Therefore, Eq. transforms as 

0{t,x)^+jiOit,-x). (2) 

This means that the intrinsic parity of Eq. ^ is pos- 
itive. Although the intrinsic parity of Eq. |^ isposi- 
tive, it couples to negative parity states as well [52]. In 
the Dirac representation, one sees that the upper (lower) 
component creates positive (negative) parity states. In 
addition, in lattice formulation with limited temporal ex- 
tension, it may create backwardly-propagating state with 
opposite parity. Since the parity projection of baryon is 
thus involved, we explain our strategy of parity projec- 
tion. 

We consider the asymptotic behavior of the zero- 
momentum projected correlator as 

G„/3(^) = ■^E<^"(*'^)^'3('^'0)>. (3) 

X 

in large t and large Nt — t region. {V and Nt refer to 
the spatial volume and the temporal size of the lattice, 
respectively.) In this region, the correlator is assumed to 
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be dominated by the lowest-lying states of each parity, 
and is decomposed into two parts in the following manner 



where the choice of "±" depends on the boundary con- 
dition of quark fields in imaginary time, i.e., for 
anti-periodic boundary condition and "— " for periodic 
boundary condition, and m_ refer to the energies of 
lowest-lying states in positive and negative parity chan- 
nels, respectively. P± = (1 ± 74)/2 work as projection 
matrices onto the "upper" and the "lower" Dirac sub- 
spaces, respectively. 

We consider the first term in Eq. |0J. In the region 
< i < Nt/2, C+e"™+* is expected to dominate the 
correlator, whereas, in the region Nt/2 <C t <C Nt, 
C_e^™"'^^*~*-' is expected to dominate the correlator. In 
the intermediate region, i.e., t Nt/2, the contributions 
from the both terms cannot be negligible. The second 
term in Eq. ^ behaves in a similar manner. In this way, 
by confining ourselves to the region ^ t <C Nt/2, and 
by applying the projection matrices P*^^^ to G{t), we can 
pick up the positive and negative parity states separately. 

One may point out that the exact parity projection 
can be achieved by summing up G^"*"' and G^~\ because 
the backwardly-propagating contributions cancel. This 
is often used to project parity for ordinary "three-quark" 
baryons at quenched level. However, in the case of penta- 
quark, this may not work even in quenched level, because 
the five-quark system can decompose into N and K so 
that N may propagate forward whereas K may propagate 
backward or vice versa [47i . 



III. NUMERICAL RESULT ON 9+ MASS 

A. Lattice parameter set 

To generate gauge field configurations, we use the stan- 
dard plaquette action on the anisotropic lattice of the size 
12^ X 96 as 

- |-— E RcTr{l-F,,(.)} (5) 
+ E ReTr{l-P,4(s)}, 



s,i<3 



where P^iy{s) £ SU(3) denotes the plaquette operator 
in the ^-i/-plane. The lattice parameter and the bare 
anisotropy parameter are fixed as /3 = 2Nc/ = 5.75 and 
7g — 3.2552, respectively, which reproduce the renor- 
malized anisotropy as ^ = a^/at = 4 48]. Adopting 
the pseudo-heat-bath algorithm, we pick up gauge field 
configurations every 500 sweeps after skipping 10,000 



sweeps for the thermalization. We use totally 504 gauge 
field configurations to construct the temporal correlators. 
The lattice spacing is determined from the static quark 
potential adopting the Sommer parameter r^^ — 395 
MeV (ro - 0.5 fm) as aj^ = 1.100(6) GeV (a^ ~ 
0.18 fm). Hence the lattice size 12^^ x 96 amounts to 
(2.15fm)'^ X (4.30fm) in the physical unit. 

For quark fields "0 and i/i, we adopt the 0(a)-improved 
Wilson (clover) action on the anisotropic lattice as 



Sf = ^tp{x)K{x,y)ip{y), 



(6) 



K{x,y) = S^^y- Kt\^{'^-li)Ui{x) 5^j^-^y 

+ {l + lA)Ul{x-l)5,_-^^y] 

+ {r + y,)ul{x-l)5^_'^^y] 

i i>j 

where Kg and denote the spatial and temporal hop- 
ping parameters, respectively. F^i, denotes the field 
strength, which is defined through the standard clover- 
leaf-type construction, r denotes the Wilson parameter. 
ce and cb denote the clover coefficients. To achieve the 
tadpole improvement, the link variables are rescaled as 
Ui{x) Ui{x)/us and U4{x) U4{x)/ut, where Ug and 
Ut denote the mean-field values of the spatial and tem- 
poral link variables, respectively (ili'H^I. This is equiv- 
alent to the redefinition of the hopping parameters as 
the tadpole-improved ones (with tilde), i.e., Kg = kg/ug 
and Kt = kt/ut- The anisotropy parameter is defined 
as 7f = kt/Hg, which coincides with the renormalized 
anisotropy ^ = a^/ar for sufficiently small quark mass 
at the tadpole- improved level jS^. For given Kg, the four 
parameters r, ce, cb and Kg/Kt should be, in principle, 
tuned so that "Lorentz symmetry" holds up to discretiza- 
tion errors of 0{a^). Here, r, ce and cb are fixed by 
adopting the tadpole improved tree-level values as 



1 



1 



CE 



CB 



(7) 



Only the value of k„/ n^- (= ^e • {ua/ur)) is tuned non- 
perturbatively by using the meson dispersion relation 
l49ll. It is convenient to define k as 



i EE ^ - 2 (7J. - 3r - 4) , 



(8) 



then the bare quark mass is expressed as mo = ^(I/k — 8) 
in the spatial lattice unit in the continuum limit. This 
K plays the role of the hopping parameter "k" in the 
isotropic formulation. For detail, see Refs. [i^lsoj . from 
which we take the lattice parameters. The values of the 
lattice parameters are summarized in Table 
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P 7G 0,3 /at a^^ [GeV] Size A''conf Ua Ut jf Kc Values of k 

5.75 3.2552 4 1.100(6) 12^ x 96 504 0.7620(2) 0.9871(0) 3.909 0.12640(5) 0.1240, 0.1230, 0.1220, 0.1210 



TABLE I: Parameters of the lattice simulation. The spatial lattice spacing is determined with r^^ = 385 MeV for the 
Sommer parameter. The mean-field values of link variables (ua and Ut) are defined in the Landau gauge. Kc denotes the 
critical value of k. 



0.1210 



0.1220 



0.1230 



0.1240 



^phys 



0.1210 0.1220 0.1230 0.1240 empirical 



1.005(2) 0.898(2) 0.784(2) 0.656(3) 0.140 

nip 1.240(3) 1.161(3) 1.085(4) 1.011(5) 0.850(7) 

niK 0.845(2) 0.785(2) 0.723(2) 0.656(3) 0.530(4) 

mjv 1.878(5) 1.744(5) 1.604(5) 1.460(6) 1.173(9) 



TABLE IL Masses of vr, p, K and N for each hopping param- 
eter K in the physical unit GeV. Kphys. — 0.1261 denotes the 
value of K which achieves m,r — 0.14 GeV. 



We adopt four values of the hopping parameter as 
K = 0.1210,0.1220,0.1230 and 0.1240, which correspond 
to niT^/nip = 0.81,0.77,0.72 and 0.65, respectively. For 
temporal direction, we impose anti-periodic boundary 
condition on all the quark fields. For spatial directions, 
unless otherwise indicated, we impose periodic boundary 
condition on all the quarks. We refer to this boundary 
condition as "standard (spatial) BC\ 

By keeping Kg ~ 0.1240 fixed for s quark, and by 
changing k — 0.1210 — 0.1240 for u and d quarks, we 
perform the chiral extrapolation to the physical quark 
mass region. In the following part of the paper, we will 
use 



(ks,k) = (0.1240,0.1220), 



(9) 



as a typical set of hopping parameters in presenting cor- 
relators and effective mass plots. For convenience, we 
summarize masses of tt, p, K and N for each hopping pa- 
rameter K together with their chirally extrapolated values 
in Table UTI Here, the chiral extrapolations of these par- 
ticles are performed with a linear function in . Unless 
otherwise indicated, we adopt jackknife prescription to 
estimate statistical errors. 

In order to enhance the low-energy contributions, we 
use a smeared source. This is achieved by employing a 
spatially extended interpolating field of the gaussian size 
p ~ 0.4 fm in the Coulomb gauge, which is obtained by 
replacing the quark fields q{x) in Eq. by the smeared 
quark fields (/smear (2^) defined as 



i-{t, x) = Afy^ exp 



\x-y\ 

2p2 



qit,y), (10) 



p-wave 
HBC 



2.723 
2.987 
2.924 



2.528 
2.809 
2.743 



2.327 
2.629 
2.558 



2.116 
2.442 
2.367 



1.440 
1.865 
1.770 



TABLE III: Numerical values of NK thresholds (Eq. lO) for 
each hopping parameter k in the physical unit GeV, when 
spatial lattice of the size is L ~ 2.15 fm. The rightmost col- 
umn indicated by "empirical" corresponds to the thresholds 
calculated for the physical masses of N and K, tum = 0.94 
GeV and tuk = 0.5 GeV. The first and the second rows show 
the NK thresholds for the s-wave and p-wave states, respec- 
tively, when the standard BC is imposed. The third row shows 
the NK thresholds for the HBC. In the HBC, the s-wave and 
p-wave NK thresholds coincide. 



B. NK scattering states for standard BC 

The NK scattering state with zero total momentum is 
obtained as J=l/2-projection of the following state: 

Nip, s)K{-p) ) ± I iV(-p, s)K{p) ) , (11) 

where "+" is for s-wave, and "— " for p-wave. s denotes 
the spin degrees of freedom of the nucleon. In finite spa- 
tial box of the size L^, the allowed momentum is quan- 
tized as 

27r 

— n, fi&J?, (12) 
L 

due to the periodic boundary condition. By neglecting 
the interaction between N and K, the energy of this state 
is approximated as 



E{p) ~ \/mj^+p^ + Jm%+p^. 



(13) 



For s-wave, the NK threshold is expressed as Eth = 
rax + 771 TV- On the other hand, for p-wave, since Eq. 111(1 
vanishes for p — 0, the NK threshold starts from ii'th = 
E{pinin) associated with the minimum momentum 



IPminI = 



(14) 



Numerical values of these NK thresholds for each hopping 
parameter k are presented in Table UTTl 



where Af is an appropriate normalization factor. For 
practical use, we extend Eq. H10() appropriately so as to fit 
a particular choice of the spatial boundary condition. In 
this paper, we present correlators with a smeared source 
and a point sink. 



C. The correlators for Q 

In Fig.n we show the effective mass plot for the posi- 
tive parity channel, which is obtained from the correlator 
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FIG. 1: The effective mass pfot of positive parity O^. Eq. @ 
is adopted as a typical set of hopping parameters. The sta- 
tistical error is obtained with jackknife error estimate. The 
solid line denotes the result of the single-exponential fit per- 
formed over the interval 25 < t < 35. The dotted line denotes 
the p-wave NK threshold energy corresponding to the spatial 
lattice size L ~ 2.15 fm. 



FIG. 2: The effective mass plot for negative parity O for the 
typical set of hopping parameters (Eq. 0). The solid line 
denotes the result of the single-exponential fit performed over 
the interval 25 < t < 35. The dotted line denotes the s-wave 
NK threshold. 



with a smeared source and a point sink adopting Eq. I^. 
The dotted hne denotes the p-wave NK threshold for the 
spatial lattice size L ~ 2.15 fm. 
The effective mass is defined as 

"«">-'°«(Gfri))- 

where G{t) denotes the correlator. At sufficiently large t, 
contributions from excited states diminish, and the cor- 
relator is dominated by a single state with energy m as 
G{t) ~ Ae~™*. Then Eq. ((T^ gives a constant effective 
mass as mefi(t) ^ m. Thus a plateau can be served as an 
indicator of the single-state saturation. One is then al- 
lowed to perform the single-exponential fit in the plateau 
region. 

In the region < t < 25 in Fig. ^ the effective mass 
nicttit) decreases monotonically, which implies that the 
higher spectral contributions are gradually reduced. We 
find a plateau in the interval 25 < i < 35, where the 
single-state dominance is expected to be achieved. Be- 
yond t 35, the data becomes too noisy. In addition, it 
is expected to receive the contribution from backwardly- 
propagating negative parity state as mentioned before. 
Hence, we simply neglect the data for t > 35, and 
perform the single-exponential fit of the correlator with 
/(<) = Ae"™* in the plateau region, 25 < t < 35. The 
solid line denotes the best-fit result. 

In Fig. 121 we show the effective mass plot for the neg- 
ative parity channel adopting Eq. @. The dotted line 
denotes s-wave NK threshold. In the small t region, i.e., 
t < 25, mcff(t) decreases monotonically, and we see a 
rather stable plateau at 25 < t < 35. The single expo- 
nential fit at 25 < i < 35 gives the solid line. 



3.5 I • • • • 1 . . . . r 




1.5 I ' ' ' ' — ' — ' — ' — ' — ' — ^ 

0.0 0.5 1.0 

[GeV^] 



FIG. 3: me for both parity states plotted against m^. The 
triangles denote positive parity, while the circles denote nega- 
tive parity. The open symbols denote the direct lattice QCD 
data, whereas the closed symbols denote the results of the chi- 
ral extrapolation. The dotted lines indicate the NK thresholds 
for p-wave(upper) and s-wave(lower) cases. 



K 0.1210 


0.1220 


0.1230 


0.1240 


/^phys. 


-1- parity 3.18(4) 
- parity 2.80(1) 


2.99(4) 
2.59(2) 


2.80(6) 
2.39(2) 


2.64(11) 
2.17(3) 


2.25(12) 
1.75(4) 



TABLE IV: The masses of 9+ for each value of k. The first 
line corresponds to the positive parity state, and the second 
line to the negative parity state. Kphys. — 0.1261 denotes the 
value of K which achieve ~ 0.14 GeV. 
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As mentioned before, k = 0.1240 is kept fixed for s 
quark, and k = 0.1210 — 0.1240 are varied for u and d 
quarks to perform the chiral extrapolation. In Fig.|3| the 
masses of positive (triangle) and negative (circle) parity 
0+ arc plotted against m^. The open symbols denote 
the direct lattice data. Since these data behave almost 
linearly in , the chiral extrapolation is performed with 
a linear function. Such a linear behavior against toJ is 
also observed for ordinary non-PS mesons and baryons 
|50| . The closed symbols denote the results of chiral ex- 
trapolation. For convenience, we show p-wave (upper) 
and s-wave (lower) NK threshold with dotted lines. 

For the positive parity state, the chiral extrapolation 
leads to me = 2.25 GeV, which is much heavier than the 
experimentally observed 0^(1540). For negative parity 
state, on the other hand, the chiral extrapolation leads 
to me = 1.75 GeV, which is rather close to the empiri- 
cal value. Our lattice QCD results thus support that the 
negative parity state is the lowest of the 0+ spectrum. 
Table summarizes masses of 8+ of both parities for 
each hopping parameter together with their chirally ex- 
trapolated values. 



IV. THE NEW METHOD TO DISTINGUISH 6+ 
AND NK 

A. The hybrid boundary condition (HBC) 

In the p-wave NK scattering states, N and K have non- 
zero minimum momenta as |pmin| = 27r/L due to the 
finiteness of the spatial box. As a consequence, the NK 
threshold is raised by a few hundred MeV depending on 
the spatial volume. In studies of positive parity we 
can utilize this volume dependence to distinguish local- 
ized resonances from NK-scattering states. In contrast, 
in studies of negative parity 6+, NK scattering states 
are of s-wave, and N and K can have zero momenta, i.e., 
|Pmin| = 0. Therefore, the NK threshold has no volume 
dependence, i.e., ttin + tuk, which is less convenient to 
distinguish 8+ from NK scattering states. 

It would be of great worth, if we could find some pre- 
scription to raise the s-wave NK threshold by changing 
the spatial volume. This can be achieved by twisting the 
spatial boundary condition in a flavor dependent manner 
as follows. We impose the anti-periodic boundary condi- 
tion on u and d quarks, whereas the periodic boundary 
condition on s quark. (See Table IVl) We will refer to this 
spatial boundary condition as ^^hybrid boundary condition 
(HBC)". 

As a consequence of the HBC, hadrons feel the bound- 
ary condition in their own ways. Since N{uud, udd) and 
K(us, ds) contain odd numbers of the u and d quarks, 
they are subject to the anti-periodic boundary condition. 
In contrast, since Q~^(uudds) contain even numbers of u 
and d quarks, it is subject to the periodic boundary con- 
dition. Recall that, due to the finite size of the lattice. 





u quark 


d quark 


s quark 


HBC 


anti-periodic 


anti-periodic 


periodic 


standard EC 


periodic 


periodic 


periodic 



TABLE V: The hybnd boundary condition (HBC) imposed on 
the quark fields. The second line shows the standard BC for 
comparison. 



the allowed spatial momenta are quantized as 

{2ni'K /L for periodic BC 

(16) 
{2ni + l)7r/i for anti-periodic BC, 

with Ui G Z. Hence, N and K acquire the non-zero min- 
imum momenta as 

IPminI = (17) 

even for the s-wave state. By comparing this with 
Eq. H14|l . we see that HBC can raise the NK threshold 
like the p-wave case. (See Eq. (dJ.) In Table UTTI we 
also show NK thresholds for the HBC. Note that, unlike 
the standard BC, s-wave threshold and p-wave one coin- 
cide in the HBC. We see, in Table Hill that s-wave NK 
threshold is raised by about 200 — 250 MeV in the range 
of 0.1210 < K < 0.1240 in our calculation. 

On the other hand, a localized resonance Q^{uudds) 
can have zero momentum as |pmin| = 0. (See Table IVTl ') 
Therefore, the shift of me comes only through the change 
in its intrinsic structure. In this case, the shift is expected 
to be less significant than the shift due to the kinematic 
reason as is the case in N and K. Now, if the spatial 
size of 8+ is sufficiently smaller than the spatial lattice 
size L, then it is safely expected that the mass of 8+ is 
insensitive to the change of the boundary condition. 

B. Numerical result on substance of 5Q state 

Before presenting the results for 8+, we have to con- 
vince ourselves that localized resonance states are really 
insensitive to the change of boundary condition. For this 
purpose, we apply the HBC to the S(uds) baryon, which 
is an established localized resonance. We take E rather 
than the nuclcon, because S contains even numbers of 
u and d quarks. Hence, E is subject to the periodic 
boundary condition, and can have zero spatial momen- 
tum. Therefore, if the mass of S is affected by the choice 
of the BC, it is attributed to a change of intrinsic struc- 
ture. In Fig. 01 we show effective mass plots of E baryon 
correlators adopting the typical set of hopping parame- 
ter Eq. The upper one corresponds to the standard 
BC, and the lower one to the HBC. We see that the 
plateau is raised by only a negligible amount due to the 
change of the boundary condition. This example explic- 
itly shows that localized resonance states are insensitive 
to the change of the boundary condition. 
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quark content 


spatial BC 


minimum momentum 






N 


uud, udd 


anti-periodic 


Pmin = (^77' "'"L ' "'"L ) 


Pmin 


= V3^ 


K 


us, ds 


anti-periodic 




|Pmin| 




e+ 


uudds 


periodic 


Pmin = (0, 0, 0) 


IPmin 


= 



TABLE VI: The consequence of the HBC on the hadrons. 



Standard BC 



> 

u 

O 



2.0 
1.9 
1.8 
1.7 



1.6 I 
1.5 



£ (uds; j''=l/2'^) 



10 



20 

tK] 

Hybrid BC 



30 



40 



> 
O 



2.0 
1.9 
1.8 
1.7 
1.6 
1.5 



L (uds; j"=l/2'^) 



10 



20 

tK] 



30 



40 



FIG. 4: The effective mass plots of Ti{uds) baryon under the 
standard BC, and under the HBC. Eq. @ is adopted as a 
typical set of the hopping parameters. The solid lines denote 
the results of the single exponential fit performed in the region 
2Q <t < 30. No significant difference is observed between the 
two results. 



Now, we present the numerical result of negative par- 
ity 8+. In Fig. |S1 we show the effective mass of nega- 
tive parity 8+ imposing the HBC adopting Kg = 0.1240 
and K = 0.1220. In the region < f < 25, the effec- 
tive mass decrease monotonically as before. There is a 
plateau in the region 25 < t < 35. Beyond t ^ 35, 
the plateau breaks down due to the contributions from 
backwardly-propagating positive parity states. Hence, 
we simply neglect the data for t < 35, and perform the 
single-exponential fit in the plateau region, 25 < t < 35. 
The solid line denotes the result of the single-exponential 
fit. T he dotted lin e denotes the modified NK threshold, 



p2 



with 



^/i^l/L 



> 

(U 

O 



CD 

s 



4.0 



3.5 



3.0 



2.5 



0+(j''=l/2") 



« Hybrid BC » 



^55 



ills 




10 



20 



30 40 



FIG. 5: The effective mass plot for the negative parity 
under the HBC. Eq. is adopted as a typical set of hopping 
parameters. The solid line denotes the result of the single- 
exponential fit performed within the interval 25 < t < 35. 
The dotted line denotes the s-wave NK threshold energy mod- 
ified due to the HBC. This figure should be compared with 
Fig.H 



due to the HBC. Comparing Fig. Elwith Fig. 12 wc ob- 
serve that the plateau is raised by about 200 MeV, which 
coincides with the shift of the NK threshold. In Fig.|Hl we 
compare the result of HBC with that of standard BC for 
each hopping parameter k for the negative parity O"*". 
Closed circles denote the results obtained with best-fit 
analysis. Solid lines show the corresponding NK thresh- 
olds. For all K, we see that the best-fit masses(~ the 
location of the plateaus) are raised by about the same 
amount as the corresponding NK thresholds. Therefore, 
we find that there is no localized resonance states be- 

= %/37r/L. 



As a result, we conclude that the negative parity state 
observed before is an NK scattering state. 



C. Comments 

(i) The fully anti-periodic boundary condition is less 
convenient for our purpose. In fact, in this case, K is 
subject to the periodic boundary condition, whereas N 
and 0+ are subject to the anti-periodic boundary condi- 
tion. Since K can have the zero momentum, the shift of 
the NK threshold is smaller. In addition, 0+ must have 
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3.0 



2.0 



- 




31 


K=0.1210 




a: 


z 


K=0.1220 
K=0.1230 
K=0.1240 


- 


a: 
Z 


I 

/ ' I 




Standard BC 


Hybrid BC 



FIG. 6: Comparison of results of starrdard BC(l.h.s.) and 
HBC(r.h.s.) for each hopping parameter k. Closed circles 
denote results obtained from the best fit analysis. The solid 
lines denote the corresponding NK threshold. 



non-zero momentum, which raises the location of 0"*" due 
to the kinematic reason. It is easy to understand that the 
resulting spectrum is less convenient for our purpose. 

(ii) It is possible to use the HBC in the study of the 
positive parity state. However, since the change of the 
minimum momenta amounts to only a minor modifica- 
tion from 2tt / L to \/iTT/L, the resulting change in the 
spectrum is less significant than the negative parity state. 

(iii) The HBC provides us with a general tool, the use 
of which is not restricted to the calculation of . For 
instance, it can be used in the study of the hadronic 
molecules. Hadronic molecules are bound state of two or 
more hadrons. The binding energies are typically 10 — 50 
MeV. In the lattice QCD studies, the gap of 10 - 50 
MeV may be too small to confirm the existence of the 
bound state. Therefore the HBC is useful by enlarging 
this gap to a few hundred MeV. 

(iv) Careful reader may point out that, besides the sta- 
ble plateau which we discussed, there is a small plateau- 
like structure in Fig. |21 and Fig. [S] in the region 15 < 
t < 18. However, this "plateau" is also raised by about 
200 MeV due to the HBC. This means that this region 
is not dominated by a single state, and a possible reason 
of this shift would be that contaminations of low-lying 
states becomes reduced in the HBC. Still, the nature of 
this small plateau cannot be determined from our data 
available so far. At any rate, we can state that it does 
not correspond to the experimentally observed 0"'"(154O), 
because its location is rather high m ~ 1.95 GeV even 
after the chiral extrapolation. Note that we do not find 
such a plateau-like structure in the effective mass plot 
with a point source and a point sink. 



V. SUMMARY AND DISCUSSION 

We have studied the penta-quark baryon 0+ in 
anisotropic lattice QCD at the quenched level. We have 
used the standard plaquette action for gauge field config- 
urations on the anisotropic lattice of the size 12'^ x 96 
with the renormalized anisotropy as /at = 4 at /3 = 
5.75, i.e., as ~ 0.18 fm, at — 0.045 fm. For quarks, 
we have adopted 0(a)-improved Wilson (clover) ac- 
tion with four values of hopping parameters as k = 
0.1210(0.0010)0.1240. To enhance the low-lying spectra, 
we have adopted a spatially extended operator. 

We have used 504 gauge field configurations to con- 
struct the correlator of a non-NK type interpolating field. 
For each parity channel, we have found a rather stable 
plateau, where the single exponential fit has been per- 
formed. After the chiral extrapolation, we have found 
me = 2.25 GeV for positive parity state, which is too 
heavy to be accepted as the experimentally observed 
0+(154O). For negative parity case, we have obtained 
TOe ~ 1-75 GeV, which is rather close to the empirical 
value. Our data have thus suggested that the negative 
parity state is the lowest of the penta-quark spectrum. 

In order to clarify whether this negative parity state 
is a localized resonance or not, we have proposed a 
new method with ^''hybrid boundary condition (HBC)" , 
where the anti-periodic spatial boundary condition is im- 
posed on the u and d quarks, while the periodic spatial 
boundary condition is imposed on the s quark. As re- 
sults, N{uud,udd) and K{us,ds), both of which con- 
tain odd numbers of u and d quarks, are subject to 
the anti-periodic spatial boundary condition, whereas 
0"*" (uudds) , which contains even numbers of the u and d 
quarks, is subject to the periodic boundary condition. A 
remarkable feature of the HBC is that it raises the s-wave 
NK threshold by a few hundred MeV without affecting 
the localized penta-quark resonance. 

By using the HBC method, we have further investi- 
gated the negative parity state. We have found that the 
plateau is raised by about 200 MeV due to the HBC, 
which is consistent with the shift of the NK threshold. 
We co nclude that there is no localized resonance state be- 
low \/m5~T^^+ y/m^ + Pmin ^it^^ \Pmin\ = VStt/L. 
It follows, in particular, that the negative parity state 
observed in the standard BC is a mere NK scattering 
state. 

Now we discuss possible implications of our results. 
Experimentally, a number of groups P, IE IE 13 have 
confirmed the existence of 0+(154O). Hence, it should be 
observed in lattice QCD Monte Carlo calculation as well. 
However, in neither positive nor negative parity states, 
we have found any relevant signals of narrow penta-quark 
states. In this way, we have arrived at a similar result as 
Ref.li^. However, it should be remarked that this result 
has been obtained in spite of our choice of the non-NK 
type interpolating field. These two results suggest that 
may not be reachable by a simple five-quark inter- 
polating field operator in lattice QCD — it might have 
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more complicated structure than expected. Needless to 
say, it is desirable to examine the finite volume artifact 
and discretization error, as usual. It is also desirable 
to reduce the ambiguity from the chiral extrapolation. 
Concerning the light quark effects, one might wonder 
that the dynamical quark effects might be important, be- 
cause the narrow resonance 0+ was originally predicted 
by the Skyrme model j^]- After all, the Skyrmion is a 
nontrivial configuration of pion. However, this argument 
is not quite correct, since the nucleon serves as a counter- 
example. A nucleon is indeed the ground state of in the 
Skyrme model and thus contains a non-trivial configu- 
ration of pion. It is however also fairly well reproduced 
in quenched lattice QCD with conventional three-quark 
interpolating field. 

Besides these rather technical stuff, the difficulty may 
trace back to our insufficient understanding of the true 
nature of the pcnta-quark. It would be interesting to 
consider the following possibilities. 

(1) Highly non-local structure: If the structure of 8+ is 
really the Jaffe-Wilczek type or the Karliner-Lipkin type, 
it may be better to invent series of new non-local inter- 
polating fields, which have sufficiently complicated struc- 
tures so that they can fit Jaffe-Wilczek/Kaliner-Lipkin 
pictures. In this respect, the group theoretical construc- 
tion of interpolating fields proposed in Ref.js^l together 
with the variational method would be interesting. 

(2) Other quantum numbers: In most lattice QCD 
calculations of 6+, spin of 6+ is assumed to be 1/2. 
Naively, one may expect that J=3/2 state would be heav- 
ier than J=l/2 states. However, since there is no exper- 
imental evidence on the spin of O"*", it would be worth 
keeping J=3/2 and 5/2 possibihties as candidates. We 
have to keep in mind that O"*" is already a quite mysteri- 
ous particle. In addition, although there is no evidence on 
the doubly charged 0++ in the experiment 8], it might 
be worth studying 1=1 and 2 states as well, which have 
been suggested by Ref.js^. Since we do not know the 
production mechanism of the penta-quark, it might be 
better to keep them as possible candidates. 

(3) KttN hepta-quark hypothesis: All the lattice QCD 
calculations for 0+ available so far are based on the 
quenched QCD. It is known that some hadrons are dif- 
ficult to be reproduced in the quenched QCD. For in- 
stance, quenched lattice QCD cannot reproduce A(1405) 
as a three-quark state [H^l- It is conjectured to be a KN 
bound state. In this case, A(1405) is a "penta-quark". 
Hence, it would be natural that mass of A(1405) can- 
not be properly reproduced in quenched lattice QCD 
with a standard three-quark interpolating field. Actu- 
ally, 0"''(154O) itself is hypothesized to be a bound state 
oiKnN, i.e., the hepta-quark [Ulllll^. If this is really 
the case, it would be natural that 0"'"(154O) is difficult 
to be observed in quenched lattice QCD with "ordinary" 
penta-quark interpolating fields. However, concerning 
the hepta-quark picture, we should keep in mind that, 
in spite of its many good features, the interaction among 
KttN system is said to be too weak to make a bound 



state with the binding energy of 40 MeV. 

In this way, it is necessary to perform more systematic 
studies of the penta-quarks in order to reveal its myste- 
rious nature. The proper use of lattice QCD can help us 
proceed in this direction. 
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APPENDIX A: REARRANGEMENTS OF EQ. dP 

The aim here is to estimate the contribution from NK- 
type operator to the correlator of Eq. To this end, 
we use the Fierz rearrangements. By using an identity 
eabcebfg = ^cfSag " ScgSaf, and by applying the Fierz 
identity to is rearranged as 

eabc^ade^bfg (u'^Cjsde) (uJCdg) Cs^ (Al) 
= Cade (uIC754) {{ulCda) Cs^ - {ulCd,) Cs^] 
= fade (u^Cjsde) 

X ^ I - J5da {Scl5Uc) + 75^0 [Scl^dc) 
-da (ScUc) + Uc (Scdc) 

+l5lf^da (Sc 757m ""c) - IZlfiUa {Sclblfidc) 
+ ^0-iiuda {ScO-^uUc) - ^(Jf^i^Ua {ScCTf^i^dc) | , 

where only the first and the second term factorize into 
the product of the interpolating fields of N and K up to 
a factor of 75 . We see that Eq. ^ contains the NK-type 
interpolating field with a factor of j. 

In Eq. (|A1|I . u and d fields in K come from the pseudo- 
scalar diquark operator. There is another contribution, 
where u and d fields in K come from the scalar di- 
quark operator. This is obtained by using the identity 
tabc^ade = Sbd^ce - Sbe^cd and by applying the Fierz iden- 
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tity to (u^Cj^d^j Cs^ in the following way: 

iabceadeebfg (u'^Cj^de) (ujcdg) Cs^ (A2) 
= e,f, (uJCrfg) { {ulC^^d,) Csl - {ulC-i,db) Csl) 

= <^bfg{uJCdg) 

X - I - u& (scjsdc) + db (scjbUc) 

-75Mb {Scdc) + l^db (ScUc) 

-JfiUb (sc 757^4) + ^i^i^db {sclbli^iUc) 



where only the first and the second terms contribute to 
the product of interpolating fields of N and K up to a 
factor of 75. The nucleon interpolating field here has an 
unfamiliar form. Another Fierz rearrangement of this 
nucleon interpolating field can make this into the famil- 
iar form, which provides an additional factor of i. To- 
tally, the contribution of the NK-type interpolating field 
to Eq. (P) in this case involves the factor of 1/16. 

By neglecting the small contribution in Eq. lA2p for 
simplicity, the direct contribution of the NK-type inter- 
polating field (in source and sink) to the correlator of 
Eq. involves the factor of 1/16, which is rather small. 
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